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ABSTRACT 

Free convection MHD flow of a viscous incompressible fluid flow through porous medium with mass and heat 
transfer past an infinite vertical plate in slip -flow regime in the presence of variable suction and periodic temperature and 
mass concentration has been studied. Assuming variable suction at the plate, approximate solutions are obtained for 
velocity, temperature, concentration, Nusselt number, Sherwood number and skin friction. The result of various material 
parameters are discussed on flow variable and presented by graphs and tables. 
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INTRODUCTION 

The study of an electrically conducting fluid, which influences many natural and manmade flows, have many 
applications in engineering problems such as magneto hydrodynamics (MHD) generators, plasma studies, nuclear reactors, 
geothermal energy extraction, and the boundary layer control in the field of aerodynamics. 

The process of heat and mass transfer in free convection flow have attracted the attention of a number of scholars 
due to their application in many branches of science and engineering, viz. in the early stages of melting adjacent to a heated 
surface, in chemical engineering processes which are classified as a mass transfer process, in a cooling device aeronautics, 
fluid fuel nuclear reactor. The phenomenon of free convection arises in the fluid when temperature and mass concentration 
change cause density variation leading to buoyancy forces acting on the fluid elements. 

The study of fluid flows is of paramount importance in engineering. For example, the unsteady free convection 
flow along a vertical plate has been given considerable interest, because of its application in devices which are cooled by 
natural convection, as in the case of electrical heaters and transformers. Joule heating in electronics and in physics refers to 
the increase in temperature of a conductor as a result of moving electrons colliding with atoms, whereupon momentum is 
transferred to the atom, increasing its kinetic energy. 

Revankar, [1] considered Natural convection effects on MHD flow past an impulsively started permeable plate. 
Mohapatra et.al [2] have been analyzed hydrodynamic free convection flow with mass transfer past a vertical plate. 
Das et al. [3] studied the transient free convection flow past an infinite vertical plate with periodic temperature variation, 
because the free convection is enhanced by superimposing oscillating temperature on the mean plate temperature. 

Hossain et al. [4] studied the influence of fluctuating surface temperature and concentration on natural convection 
flow from a vertical flat plate. Free convective flow through a porous medium between two vertical parallel plates has been 
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studied by Singh [5]. The flow regime is called the slip flow regime and this effect can not be neglected. Using this effect 
Sharma and Chaudhary [6] studied the effect of variable suction on transient free convective viscous incompressible flow 
past a vertical plate with periodic temperature variations in slip-flow regime. Recently, Sharma and Sharma [7] studied 
influence of variable suction on unsteady free convective flow from a vertical flat plate and heat transfer in slip-flow 
regime. Sharma [8]have studied Influence of periodic temperature and concentration on unsteady free convective viscous 
incompressible flow and heat transfer past a vertical plate in slip-flow regime. 

It is proposed to study the Magnetic effects of mass and heat transfer on free convection flow through porous 
medium past an infinite vertical plate in slip -flow regime in the presence of chemical reaction, variable suction and 
periodic temperature and mass concentration 

FORMULATION OF THE PROBLEM 

An unsteady free convection two dimensional flow of an incompressible and electrically conducting viscous fluid 
through porous medium past an infinite vertical plate in slip-flow regime in the presence of chemical reaction and variable 
suction is considered. The X -axis is taken along the plate in vertical upward direction and Y -axis is taken normal to it. 
A magnetic field of uniform strength H is applied normal to the plate. Initially surrounding fluid are at rest and are the 

temperature T x and mass concentration at all points. Also the plate temperature and mass concentration are 
respectively T w and C w . Since the plate is considered infinite along X direction, all physical quantities will be independent 
ofx . Under these assumptions, the physical variables are function of y and t only. Then neglecting viscous dissipation and 
assuming variation of density in the body force term. Then by usual boussinesq's approximation the unsteady flow 
governed by following equations 

Equation of continuity 

which is satisfied with v — — v n variable suction / injection 

The momentum equation 

% - v (l + eAe^') | = v^ + gP 0" - O + gfi e (£' - O - f - ^ (2) 
The energy equation 

-T- v„ 1 + eAe l01t )—, = j (3) 

at °V J dy pC pdy ' 2 v ' 

The mass concentration equation 

£- Vo ( 1 + eAe^y£ = D d ^-RXc'-C) (4) 

Where p is the density, vis the kinematic viscosity, g is the acceleration due to gravity, ^4is the suction / injection 
parameter, /?is the coefficient of volume expansion for heat transfer, /? c is the coefficient of volume expansion for mass 
transfer, D is the chemical molecular diffusivity, k is the thermal conductivity, K is the permeability of porous medium, 
a is the electrical conductivity, C p is the specific heat at constant pressure,/?' is chemical reaction parameter, B = H [i e is 
and the other symbols have their usual meanings. 
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With corresponding boundary conditions 

y'=0; u' = L' d ^,,T' = T W + e(T w -T^e^'.C' = C w + e(C w - Oe toV 
Asy'^ 0:u'^ 0,T'-»tC'-* C 
We now introduce the following non-dimensional quantities into equations (2) to (5) 
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(5) 



' V 4v' v ' T w —T^ C w -C, 



Gr = """■ T r T "\ em = •"">'<- c -- c -\ K = = °-^,h = <a. = 



4va> 



(6) 



where Pr is the Prandtl number, 5c the Schmidt number, Gr the Grashof number for heat transfer, Gm the 
modified Grashof number for mass transfer, R the chemical reaction parameter, K the permeability parameter porous 
medium, h rarefaction parameter and M the Hartmann number. 



With the help of equation (6) the equations (2) to (5) reduce to 
+ eAe u^ _ = _ +Grd+grn(p ___ Mu 

4 dt v J dy Prdy 2 

4 dt v ' dy Sc dy 2 ^ 

and boundary conditions to the problem in the dimensionless form are 

y = O-.u = h d ^,6 = 1 + ee iM ,<p = 1 + ee iM 
y -> oo : u -> 0, -> 0, <p -> 

SOLUTION OF THE PROBLEM 



(7) 
(8) 
(9) 



(10) 



Assuming small amplitude oscillations (e« 1), we can represent the velocity u, temperature 0and concentration 
(p near the plate as follows 



u = u (y) + eu 1 (y)e 1<0£ 
= o (y) + e0 1 (y)e i<at 
<P = <Po(y) + e(p r {y)e i<Dt 



(11) 



Substituting (11) in (7) to (10), equating the coefficients of harmonic and non harmonic terms, neglecting the 
coefficients of e 2 , we get 

u'q + u - Q; + m) u = -Gr8 - Gm<p x 
+ u[ - + ^ + M ) u ± = -Au - Gr8 1 - Gm(p 1 



8 + Pr8 = 

8'; + Pr8' 1 - P -^=-APr8 

(p + Sc<p — ScRq> = 
<pl + Sc(p[ - + RSc^ <p 1 — -ASc(p 



} 



(12) 
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With the following boundary conditions 
u = h d -^,u ± = h d -^,9 = = l,(p = l ><Pl = laty = 



, Pr+yjpr 2 +Pria> Sc+V Sc 2 +4ScR Sc-i-J Sc 2 + u>iSc+4ScR 

where a t = , a 2 = , a 3 = , 



0-4 = ~ a 5 



A 2 = ; ~ ASC l 2 isc „ Mi = 1 - A 2 A 3 = 



a 2 -Sca 2 -^-RSc' 1 2 ' 3 Pr2-Pr-(i + M) 

_ -Gm _ -(hA 3 Pr+hA 4 a 2 +A 3 +A 4 ) 

A * ~ a2 2_ a2 _(± +M y A 5 ~ ^ 



4PrGrA „ /„ 4Prjl\ 
/L4,Pr : -Gr 1+ — 

f 10) /l V 0)1 ) 

p r z_p r _(^ + I +M )< 7 a z_ ai _ ( ^ + i +M) ' 



. _ AA 4 a 2 -GmA 2 . -GmA 



_ ^ 5 a 4 
"10 — 



Al = 



«i-«4-(f4 +M ) 

-(h/l 6 Pr+ha 1 /l 7 + ha2/l8 + ' Ia 3'49 + ' la 4' 4 io+'46+'l7+'l8+'49+'4io) 



(13) 



u = 0, Ui = 0, O = 0, d 1 — 0, (p = 0, (p ± — as y -» oo 
Solving the set of equation (12) under the boundary conditions (13) we get 

O = e-T (14) 

X = fl + — ^e-^-— e" Pry (15) 

<Po = e-^ (16) 

<p x =A 1 e~ a3y +A 2 e~ a2y (17) 

u = A 5 e- a * y + A^e'^y + A 4 e~ a2y (18) 

Ul = A lie -^y + A 6 e- pi ~y + A 7 e~ a iy + A 8 e- a *y + A 9 e~ a iy + A 10 e- a *y (19) 



Let M r and M t , 9 r and 6 t and <p r and (p t are the real and imaginary parts of u 1 9 1 and <p 1 respectively. 
Now equation (19) can be written as u ± (y, t) = M r + iM t 
Also u(y, t) = u (y, t) + eu r (y, t)e ia>t 

= [u + e(M r cosa>t — Misincot)] + ie[M r sina>t + Micoscot] (20) 
The real part of u(y, t) 

u r — u + e(M r cosa>t — Misincot) (21) 
This expression takes simplest form for a)t — | and given by u r = u — eM t (22) 
Again equations (15) and (17) can be written as 0i(y, t) = r + id t and C ± (y, t) = (p r + iq>i 
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The real part of 9(y, t) 

6 r = 6> + e[6 r sina>t + Oicoscot] (23) 
This expression takes simplest form for a)t — - and given by 6 r — 9 — eQi (24) 
The real part of <p(y, t) 

<p r — cp + e[(p r sina>t + (picoscot] (25) 
This expression takes simplest form for a>t — - and given by (p r — <p — ecpi (26) 
In view of equations (1 1), (18) and (19), non-dimensional the skin friction 

T °-ye y ) y=0 -{ ay +ee 3y) y=0 

= -(a 4 ^ s + PrA 3 + a 2 A 4 ) + ee i0Jt (F r + iF t ) 
= (t to ) + eFe i(w£+a) 

= T m + e|F|e i(wt+a) (27) 

Where \F\ — J F? + Ff , tana — — , x m —— -(a 4 A 5 + PrA 3 + a 2 A 4 ) is skin friction of mean velocity and 
F — F r + iFj = -(a 5 A 11 + PrA 6 + a ± A 7 + a 2 A 8 + a 3 A g + a 4 v4 10 ) 
In view of equations (1 1), (14) and (15), the Nusselt number 

Nu = -(ir) 

= (Pr)-ee itat (S 1 + iS 2 ) (28) 
where \Q\ = Jsf+S*, tan(3 = S f and Q = S ± + iS 2 = -a r (l + ^fj + 
In view of equations (11), (16) and (17), the Sherwood number 

sh = -m 



v d y d y ' y =o 

= (a 2 ) - ee iMt (S 3 + iS 4 ) (29) 



where |r| = JS 3 + 5|, tany = — and T = -(a 3 A ± + a 2 A 2 ) 

RESULTS AND DISCUSSIONS 

In this paper we have studied the Magnetic effects of mass and heat transfer on free convection flow through 
porous medium past an infinite vertical plate in slip -flow regime in the presence of chemical reaction, variable suction and 
periodic temperature and mass concentration. The effect of Gr, Gm, M, K, R, Pr, Sc, A, h, <w and t on flow characteristics 
have been studied and shown by means of graphs and tables. 
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Velocity Profiles 

The velocity profiles are depicted in Figures 1-4. Figure 1 illustrates the effect of the parameters R, K and h on 
velocity at any point of the fluid, when Gr=l, Gm=l, K=l, Sc=l, Pr=l, t=l, A=l and co=l. 

It is noticed that the velocity increases with the increase of rarefaction parameter (h) and porous parameter (K), 
where as decreases with the increase of chemical reaction parameter (R). 

Figure 2 illustrates the effect of the parameters Gr, Gm and M on velocity at any point of the fluid, when Sc=l, 
Pr=l, R=3, K=l, t=l, h=l, A=l and 0)=1. It is noticed that the velocity decreases with the increase of magnetic parameter 
(M) and modified Grashof number (Gm) where as the velocity increases with the increase of Grashof number (Gr). 

Figure 3 illustrates the effect of the parameters Sc, t, a> and Pr on velocity at any point of the fluid, 
when Gr=l, Gm=l, R=3, K=l, M=l, h=l and e=0.2. It is noticed that the velocity decreases with the increasing of Schmidt 
number (Sc), prandtl number (Pr) and oscillating frequency (a>), where as decreases with the increase of time (t). 

Figure 4 illustrates the effect of the parameter A on velocity at any point of the fluid, when 
Gr =1, Gm=l, R=3, K=l, M=l, h=l, Sc=0.3, Pr=0.71, t=l, a> = 1 and e=0.2. It is noticed that the velocity increases with 
the increase of suction parameter(A). 

Temperature Profile 

The temperature profiles are depicted in Figure 5 only. Figure 5 illustrates the effect of the parameter Pr, A, 
t and a) on temperature at any point of the fluid at e — 0.2, in the absence of other parameters. It is noticed that the 
temperature falls with of increasing of prandtl number (Pr) and suction parameter (A), where as it falls initially then rises 
with the increase of time (t) and oscillating frequency (oj). 

Mass Concentration Profile 

The mass concentration profiles are depicted in Figure 6 only. Figure 6 illustrates the effect of the parameters Sc, 
A and R on mass concentration at any point of the fluid at a> — 1, t=l and e = 0.2, in the absence of other parameters. 

It is noticed that the mass concentration decreases with the increasing of Schmidt number (Sc), reaction parameter 
(R) and suction parameter (A). 

Shearing Stress/Skin Friction 

Table 1 illustrates the effect of the parameters A, Pr, Sc, Gr, Gm, R, M, K, h, t and co on Absolute value of skin 
friction(x), skin friction of mean velocity (t ), Amplitude of skin friction(F) and Phase angle tan(a). It is noticed that 

• The phase tan (<p) of skin friction at plate increases with the increase of magnetic parameter (M), modified grash 
of number (Gm), oscillating frequency (a>) and chemical reaction parameter (R), and decreases with the increase 
of rarefaction parameter (h). 

• The Absolute value of skin friction(x) at plate increases with the increase of suction parameter (A), porous 
parameter (K), grashof number (Gr), modified grashof number (Gm), prandtl number (Pr) and Schmidt number 
(Sc), and decreases with the increase of chemical reaction parameter (R), magnetic parameter (M), rarefaction 
parameter (h) and oscillating frequency (o>). 

• The skin friction of mean velocity (t ) at plate increases with the increase of porous parameter (K), grash of 



Magnetic Effects of Mass and Heat Transfer on Free Convection Flow through Porous Medium 25 
Past an Infinite Vertical Plate in Slip - Flow Regime in the Presence of Chemical Reaction, 
Variable Suction and Periodic Temperature and Mass Concentration 

number (Gr) and modified grashof number (Gm), and decreases with the increase of rarefaction parameter (h), 

chemical reaction parameter (R), prandtl number (Pr), Schmidt number (Sc) and magnetic parameter (M). 

• The Amplitude of skin friction(F) at plate increases with the increase of suction parameter (A), porous parameter 
(K), grashof number (Gr), modified grashof number (Gm), prandtl number (Pr) and Schmidt number (Sc), 
and decreases with the increase of magnetic parameter (M), oscillating frequency (a>) and rarefaction parameter 
(h). 

Nusselt Number 

Table 2 illustrates the effects of the parameters A, Pr, a) and t on Absolute value of Nusselt Number(Nu a ), 
Amplitude of Nusselt Number(Q) and Phase angle tan(/?). It is noticed that 

• The Absolute value of Nusselt Number(Nu a ) at plate increases with the increase of suction parameter (A) and 
prandtl number (Pr), and decreases with the increase of oscillating frequency (a>) and time (t). 

• The Amplitude of Nusselt Number(Q) at plate increases with the increase of suction parameter (A) and prandtl 
number (Pr), and decreases with the increase of oscillating frequency (a>). 

• The Phase angle tan( if/ ) at plate increases with the increase of oscillating frequency (<d) and prandtl number (Pr), 
and decreases with the increase of suction parameter (A). 

Shreewood Number 

Table 3 illustrates the effects of the parameters A, R Sc, a> and t on Absolute value of Sherwood Number (Sh a ), 
Amplitude of Sherwood Number (T) and Phase angle tan(y). It is noticed that 

• The Absolute value of Sherwood Number (Sh a ) at plate increases with the increase of suction parameter (A), 
Schmidt number (Sc), chemical reaction parameter (R), and decreases with the increase of oscillating frequency 
(gj) and time(t). 

• The Amplitude of Sherwood Number (T) at plate increases with the increase of suction parameter (A), Schmidt 
number (Sc) and chemical reaction parameter (R), and decreases with the increase of oscillating frequency (oj). 

• The Phase angle tan(y) at plate increases with the increase of, chemical reaction parameter (R) and oscillating 
frequency (a>), and decreases with the increase of suction parameter (A) and Schmidt number (Sc). 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.5 1 .8 2 

Figure 1: Effect of R, K and h on Velocity Profile When 
Sc=l, Pr=l, Gr=l, Gm=l, t=l, M=l and o=l 
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Figure 2: Effect of Gr, Gm and M on Velocity Profile When 
Sc=l, Pr=l, R=3, K=l, t=l, h=l and <w=l 




Figure 3: Effect of Sc, Pr <w and t on Velocity Profile When 
Gr=l, Gm=l, R=3, K=l, M=l, h=l and e=0.2 




Figure 4: Effect of A on Velocity Profile When Gr=l, Gm=l, 
R=3, K=l, M=l, h=l, Sc=0.3, Pr=0.71, t=l, o = 1 and e=0.2 
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□ 0.2 0.4 O.C □.□ 1 1.2 1.4 1.G 1.0 2 

Figure 5: Effect of A, Pr, t and <w on Temperature Profile When e=0.2 




Figure 6: Effect of A, Sc and R on Mass Concentration Profile When e=0.2, t=l and <w = 1 



Table 1: Effect of A, Pr, Sc, Gr, Gm, R, M, K, h, t and <w on Absolute Value of Skin Friction (t), 
Skin Friction of Mean Velocity (t ), Amplitude of Skin Friction(F) and Phase Angle Tan(a) 



Initially Value of the Parameters 


Skin Friction 


Absolute 


Amplitude 


Phase 


A=2, R=3, M=l 


K=l, Gr=l, Gm=l, 


of Mean 


Value of Skin 


of Skin 


Angle 


Pr=0.71,Sc=C 


.3 )<y = l,Ec=0.2 


Velocity (t„) 


Friction(r a ) 


Friction(F) 


Tan(a) 




2 


0.3529 


0.4885 


0.7569 


-44.7983 


A 


5 


0.3529 


0.6801 


1.8525 


7.9732 




7 


0.3529 


0.8164 


2.5895 


6.4706 


R 


5 


0.3348 


0.4598 


0.6498 


-3.9838 


7 


0.3228 


0.3916 


1.1616 


0.5938 


M 


3 


0.2285 


0.3232 


0.4837 


-3.0772 


5 


0.1726 


0.2774 


0.3631 


-2.1768 


K 


3 


0.4409 


0.5927 


0.9704 


4.2797 


5 


0.4653 


0.6185 


1.0368 


3.1425 


h 


2 


0.2117 


0.2929 


0.4540 


-59.5677 


4 


0.1176 


0.1626 


0.2522 


-76.3329 


Gr 


3 


0.7427 


1.6695 


5.0940 


11.4665 


5 


1.1326 


2.8643 


9.4359 


10.4130 


Gm 


3 


0.6688 


0.4880 


2.1105 


4.0146 


5 


0.9846 


0.7463 


4.9553 


4.8363 


Pr 


1 


0.3246 


0.5380 


1.2066 


14.9599 


7 


0.1996 


0.8268 


3.2495 


-61.2345 


Sc 


0.6 


0.3196 


0.3524 


0.9213 


0.5078 


0.78 


0.372 


0.8392 


4.8722 


-0.0329 


0) 


1.5 


0.3529 


0.4485 


0.4780 


-14.6083 


2 


0.3529 


0.4120 


0.3392 


-7.6540 
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Table 2: Effect of A, Pr, t and <w on Absolute Value of Nusselt Number (Nu a ), 
Amplitude of Nusselt Number(Q) and Phase Angle Tan(/?) 











Absolute Value 


Amplitude of 


Phase 




Pr 


a> 


t 


of Nusselt 


Nusselt 


Angle 












2 


0.71 


1 




0.9984 


1.9894 


-0.0625 


5 


0.71 


1 




1.3550 


3.3616 


-0.1658 


7 


0.71 


1 




1.5984 


5.1138 


-0.1930 


2 


1 


1 




1.4167 


2.8746 


-0.0583 


2 


7 


1 




9.4457 


20.9749 


-0.0118 


2 


0.71 


1.5 




0.8475 


1.9156 


-0.0533 


2 


0.71 


2 




0.6624 


1.8625 


-0.0331 


2 


0.71 


1 


1.5 


0.8588 


1.9894 


-0.0625 


2 


0.71 


1 


2 


0.6781 


1.9894 


-0.0625 



Table 3: Effect of CO >Sc, and t on Absolute Value of Sherwood Number(Sh a ), 
Amplitude of Shreewood Number (T) and Phase Angle Tan (a) 



A 


R 


Sc 




t 


Absolute Value 


Amplitude of 


Phase 


2 


3 


0.3 


1 




1.2869 


1.4582 


0.0219 


5 


3 


0.3 


1 




1.3624 


1.9778 


0.0108 


7 


3 


0.3 


1 




1.4148 


2.3242 


0.0062 


2 


5 


0.3 


1 




1.5925 


1.7208 


0.0153 


2 


7 


0.3 


1 




1.8420 


1.9381 


0.0118 


2 


3 


0.6 


1 




1.9710 


2.4065 


0.0166 


2 


3 


0.78 


1 




2.3375 


2.9421 


0.0145 


2 


3 


0.3 


1.5 




1.1588 


1.4592 


0.0329 


2 


3 


0.3 


2 




1.0114 


1.4607 


0.0438 


2 


3 


0.3 


1 


1.5 


1.1618 


1.4582 


0.0219 


2 


3 


0.3 


1 


2 


1.0177 


1.4582 


0.0219 
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